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ABSTRACT 

We present the jump conditions for ionization fronts with obHque magnetic fields. 
The standard nomenclature of R- and D-type fronts can still be applied, but in the 
case of oblique magnetization there are fronts of each type about each of the fast- and 
slow-mode speeds. As an ionization front slows, it will drive first a fast- and then a slow- 
mode shock into the surrounding medium. Even for rather weak upstream magnetic 
fields, the effect of magnetization on ionization front evolution can be important. 
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1 INTRODUCTION 

The sequence of evolution of the Hii region around a star 
with a strong ionising radiation field which turns on rapidly 



is well known (Kahn 1954; Goldsworthy 1961; see also 
Spitzer 1978, Osterbrock 1989 or Uyson & Wilfiams 1997). 
Initially, the ionization front (IF) between ionized and neu- 
tral gas moves outwards at a speed limited by the supply of 
ionizing photons, but it begins to decelerate as the ionizing 
flux at the front surface is cut by geometrical divergence and 
absorption by recombined atoms within the front. Eventu- 
ally, when the speed of the front decreases to roughly twice 
the sound speed in the (now highly overpressured) ionized 
gas, a shock is driven forwards into the neutral gas ahead 
of the front. Before this stage, the front is referred to as 
R-type, while subsequently it is referred to as D-type. The 
shock propagates outwards, gradually weakening, until, in 
principle, the Hn region eventually reaches pressure equi- 
librium with its surroundings. 

Where the external medium has an ordered magnetic 
field, the obvious critical flow speeds are the fast, Alfven and 
slow speeds rather than the isotherm al sound speed. Red- 
man et al. (1998, hereafter Paper I) studied IFs with the 



magnetic field vector in the plane of the front, and found 
that the fast-mode speed plays the same role as the sound 
speed does in the hydrodynamic case. In this paper, we ex- 
tend their work to treat the case of an IF moving into a 
medium in which the magnetic field is oblique to the di- 
rection of propagation of the front (note also we here follow 
the more conventional usage in which the magnetic fields are 
termed parallel or perpendicular with respect to the front- 
normal). 

Jump conditions for IFs with oblique magnetization 



have previously been studied by Lasker ( 1966 ) . Here we con- 
sider a wider range of upstream magnetic fields, since obser- 
vations have shown that higher magnetic fields a re found 
around H n regions than once thought likely (e.g. |R.obert 



et al. 199^; Roberts et al. 1995). We determine the properties 



of the jumps as functions of upstream conditions, using the 
velocity of the front as a parameter rather than as the vari- 
able for which we solve. We use evolutionary conditions to 
isolate the stable IF solutions, and verify these conclusions 
for a simple model of the internal structure of the fronts 
and using numerical simulations. We find that rather weak 
parallel magnetic fields can lead to a substantial decrease 
in the D-critical (i.e. photoevaporation) velocity from dense 
clumps except where the magnetic field is exactly parallel 
to the IF, and also find additional solutions to the jump 
conditions in the range of front velocities forbidden by the 
hydrodynamical jump conditions, which were not considered 
by Lasker. 

In the following sections, we present the jump condi- 
tions for MHD shocks (Section and discuss the regions 
for which evolutionary conditions suggest that these solu- 
tions are stable (Section ^). We verify that the evolutionary 
solutions are those with resolved internal structures for one 
simple model for the internal structure of the fronts (Sec- 
tion ^). In the context of these results, we discuss the de- 
velopment of an IF over time (Section p|) and illustrate this 
development using numerical models (Section P). 

Finally (in Section |^, we summarize our results, and 
provide an example of their application to observations of 
the Hn region S106. Our physical interpretation of the de- 
velopment of MHD IFs will, we hope, facilitate the future 
application of these results. 
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2 JUMP CONDITIONS 

We orient axes so that z is normal to the front, and (without 
loss of generality) that the upstream velocity and magnetic 
field are in the {x,z) plane. We use the usual MHD jump 
conditions: 



[pv,] = 



-0- 



pvzVx J- D^DxlAi\] — ^ — e- 
[Bz] = 

[VxB:,—VzBx\ = 0, 



(1) 



(4) 
(5) 



except that instead of using the energy flux condition, we 
adopt the isothermal equation of state p = pCg where the 
sound speed, Cs, increases across the front but is constant 
on either side of it. We use subscripts 1 and 2 to denote up- 
stream and downstream parameters, respectively, and write 
Vx = Ml, 2, Vz = wi,2, Cs = ci,2 and Bx = -Bi,2. Hence 



5 (6) 

p^(vl + c\) + Bl/i^ (7) 



Pl/P2 ^ V2/V1 
P2(v2+cl)+Bl/S-K 

P2W2M2 — BzB2/47r = piviui — B^Bi/Atv 

U2Bz — V2B2 = uiBz — viBi 



(8) 
i9L 



Equatii na (|[), (|) and ^ give 



B2 = 



ml 



2771 



Smf — 2r]i 



Bi, 



(10) 



where m — v/cs and we deflne r) — B^/8iTpc^ and ^ = 
B^/iiTpc? (the z and x contributions to the reciprocal of the 
plasma beta). The dependence on the upstream transverse 
velocity has disappeared, as expected as a result of frame- 
invariance. 

In equation (Q), we substitute with ( p^ for B2 and use 
equation (^ to eliminate p2 and 112 to find that, so long as 
(5 7^ and 5m{ 7^ 2r;i, the dilution fac tor 5 is given by the 
quartic equation (see also Lasker 1966) 

m?5*-m^(l + m?+4?7i+Ci)<5' 
+m\{am\ + 4r7i(l + ml + -qx + 



+ - 4771(771 + 7711(771 + Ci + ct)))5 



0, 



(11) 



+4a77i 

where a = (c2/ci)^ (100 is a typical value). It is easily ver- 
ified that this equation has the correct form in the obvious 
limiting cases (of isothermal MHD shocks where a = 1, and 
perpendicular- magnetized IFs where 771 = 0, see Paper I). 



3 SOLUTIONS 

Equation ( pj| ) relates the dilution factor, S, to the upstream 
properties of the fiow. Its roots can be found by standard 
techniques. If the flow is to have a unique solution based 
only on initial and boundary conditions, then some of these 
roots must be excluded. It is possible to exclude roots on the 
basis that they correspond to flows that do not satisfy an 
evolutionary condition, analogous to that long used in the 
study of MHD shocks. The evolutionary condition is based 
on the requirement that the number of unknowns in the 
jump conditions matches the number of constraints applied 
to the flow (e.g. [jefi'rey fc Taniuti 196^ . For MHD shocks. 



the number of characteristics entering the shock must be 
two greater than the number leaving it (since the number 
of independent shock equations is equal to the number of 
conserved variables and if there are no internal constraints 
on the front structure). 

The applicability of the evolutionary conditions to 
MHD shocks has been the subject of controversy in the 



recent pa st (e.g. Brio & Wu 1988; Kennel, Blandford & 



Wu 1990), with various authors suggesting that interme- 



diate shocks (i.e. shocks which take the flow from super- to 
sub- Alfven speeds) may be stable. However, Falle & Komis- 
sarov (1997; 1999) have shown that the non-evolutionary 



solutions are only stable when the symmetry is artiflcially 
constrained, so that the magnetic flelds ahead of and behind 
the shocks are precisely coplanar. In any cases in which the 
boundary conditions differed from this special symmetry, 
the solutions including non-evolutionary shocks were found 
to be unstable. 

The equations governing the dynamics across an IF are 
no longer a system of hyperbolic conservation laws, since the 
ionization source term cannot be neglected on the scale of 
the front. It seems reasonable to apply analogous evolution- 
ary and uniqueness conditions, but the mathematic al proofs 
for hyperbolic co nservation laws with dissipation (Falle & 



Komissarov 1999) no longer apply. The 'strong' evolution- 



ary conditions suggest that for IF the number of incoming 
characteristics is equal to the number of outgoing charac- 
teristics, since applying the ionization equation means that 
there is an additional constraint on the velocity of the front. 

Solutions which are under-specified by the external 
characteristic constraints, termed 'weakly evolutionary', 
may occur if there are internal constraints, as is the case for 
strong D-type IFs in hydrodynamics. Stable weakly evolu- 
tionary solutions only occur for limited classes of upstream 
parameters which depend on the internal structure of the 
fronts. Where the number of characteristics entering the 
front is greater than suggested by the evolutionary condi- 
tions, the solution can be realized as an MHD shock leading 
or trailing an evolutionary IF. 

In the limit in which a tends to unity from above, 
the phase change through the IF has no dynamical con- 
sequences, and the IF jump conditions approach those for 
isothermal MHD shocks. This can be seen if we rewrite equa- 
tion (0) as 

{5 - l)[m^5'^ ~ 771^(1 -f 477 + 
W(477(l + 7? + - - 477'] 

= -(a - l)[m^5 - 2??]^ 



(12) 



The IF solution which satisfies the strong evolutionary con- 
ditions becomes the trivial [5 — 1) solution of the MHD 
jump condition. The other solutions to the IF jump con- 
ditions become non-evolutionary or evolutionary isoth erma l 
MHD shocks. The analysis of Falle & Komissarov (199t) 
rules out the former as physical solutions. The latter are 
treatable as separate discontinuities, which will propagate 
away from the IF when the fiow is perturbed, since the co- 
incidence between the speeds of the shock and of the IF will 
be broken. This argument by continuity supports the use of 
the evolutionary conditions for IFs. 

As a result of this discussion, we will proceed for the 
present to isolate solutions in which the number of charac- 



© 0000 RAS, MNRAS 000, 



Magnetic ionization fronts II 3 




Figure 1. Allowed obliquely magnetized IFs (shown grey). The solid lines show the edges of the forbidden region for r; = 0, and the 
dashed lines where the inflowing gas moves at the slow, Alfven and fast speeds, from left to right. The grey regions correspond to allowed 
weak solutions for (a) r}i = 0.01, (b) rji = 3xl0~^, (c) Jji = 1, and (d) rji = 1000. In (a) the regions are labelled: 1 - the 1^1, fast-R 
solutions; 2 - the 2 —> 2, fast-D solutions; 3 — the 3 — > 3, slow-R solutions; 4 — the 4 — > 4, slow-D solutions. The fast-D and slow-R 
regions are separated by the dashed line at the upstream Alfvc5n speed. Critical solutions will be expected at the right-most edge of the 
slow-D type region (a 4 — > 3 IF) and of the fast-D type region (a 2 — > 1 IF). 



teristics entering an IF is equal to the number leaving it (and 
discuss in Section ^ the weakly evolutionary solutions for a 
simplified model of the internal structure of IFs) . The veloc- 
ity of fronts obeying the strong evolutionary conditions must 
be between the same critical speeds in the upstream and 
downstream gas (somewhat confusingly, the fronts which 
obey the strong evolutionary conditions are termed 'weak' in 
the standard nomenclature of detonations and IFs). We fol- 
low the usual classification of fiow speeds relative to the fast, 
Alfven and slow mode speeds 1 > Wf > 2 > Wa > 3 > > 4, 
so the allowed fronts are 1^1, 2^2, 3— *-3 and 4-^4. 
By analogy with the nomenclature of non-magnetized IF, 
we call these fast-R, fast-D, sIow-R and sIow-D type IF, re- 
spectively. 

The panels of Fig. |^ show regions of mi and space 
corresponding to evolutionary MHD IFs for several values 
of rji . In the figures, we see regions corresponding to the two 
distinct classes of R- and D-type solutions. The flows into 
the R-type fronts are super-fast or super-slow, while those 
into the D-type fronts are sub-fast or sub-slow. At the edges 
of the regions of solutions either the velocity into the front is 
the Alfven speed or the exit velocity from the front is equal 
to a characteristic speed (i.e. the fast mode speed at the 
edge of the fast-R region, etc.). For comparison, the solid 



lines on these plots show the edges of the forbidden region 
for perpendicular magnetization (771 = 0, see Paper I): these 
lines reach ^1 = at the edges of the forbidden region for 
unmagnetized IF, 0.05 mi ^ 20. 

Since equation ( pT| ) is a cubic in mf, quadratic in r/i and 
linear in a and ^1, there is no simple analytic form for the 
boundaries of the regions. However, certain critical values 
can be determined analytically. For 771 ^ 2a, the slow-R- 
critical line terminates where it hits the Alfven speed at 
^1 = a — 1, while the position at which the fast D-critical 
line terminates is given by 



(l + 2r,i+ei) 



■ 8Q771 



(13) 



These points are linked, respectively, by steady switch-off 
and switch-on shocks to the points on the limiting slow-D 
and fast-R critical loci at which these loci hit the axis ^1 = 0. 
The intercept between the slow-D critical locus and the axis 
is at 



277, 



Q — 2771 ' 



(14) 



for 2771 < (a — — «), beyond which the limiting value 
is mf ~ l/(4a) as for D-critical hydrodynamic IF. To il- 
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lustrate the reason for this change in solution, we rewrite 
equation (pj|) for = as 



{mis - 2r)f{ml5'^ - (1 + ml)5 + a) = 0. 



(15) 



The flow leaving a front with no upstream perpendicular 
component of magnetic field can be either at the Alfven 
speed or at the velocity of the corresponding non-magnetized 
IF. Where the flow is in the region beyond the edge of the 
slow-D-critical region shown in Figure hi (a) or (b) , the root 
of equation ( |l5| ) for flow out at the Alfven speed is not a 
real solution, since satisfying equation (^) would require that 
B| < 0. 

Even for as small a ratio between magnetic and ther- 
mal energy upstream of the front as implied by r]\ = 0.01, 
the effect of parallel magnetization on the slow-D-critical 
velocity is dramatic. Only once 771 1/(8q:) (so the Alfven 
speed in the upstream gas is below the unmagnetized D- 
critical speed) does the fast-critical locus reach ^1 = 0, so 
that the parallel magnetization may be ignored. As 771 in- 
creases, the vertical line at the Alfven speed moves across 
the plot (see Figure decreasing the region of fast-mode 
IFs and increasing that of slow- mode IFs. When 771 ^ cxd 
(a very strong parallel magnetic fleld), the (slow- mode) for- 
bidden region is identical to that in the unmagnetized case 
(independent of ^i). 

If the upstream flow is at the Alfven velocity, m\ = 
2771 , then the physical solution (when due care is taken with 
the singularity of equation (|lo|)) is often at 5 = 1, i.e. the 
ionization of the gas does not change the flow density and 
it remains at the Alfven speed. Both classes of D-type front 
have ^ > 1 (rarefy the gas), while both classes of R-type 
have 5 < 1 (compress it). The perpendicular component of 
the magnetic field, Bx, increases in a fast-R- or slow-D-type 
IF, while it decreases in a fast-D- or slow- R-type. 

We will now study the internal structure of the fronts 
for one simple model. 



4 RESOLVED FRONTS 

Up to now we have assumed, by investigating the jump con- 
ditions, that the processes within the IF take place on scales 
far smaller than those of interest for the global flow prob- 
lem. In fact, the flow structure within an ionization front 
will vary smoothly on scales comparable to the ionization 
distance in the neutral gas. The flow may take several re- 
combination lengths behind the front to relax to equilibrium. 
Here we describe the internal structure of MHD IF in one 
simple approximation, that the temperature of the gas varies 
smoothly through the front but the flow obeys the inviscid 
MHD equations througho ut (as used in the study of hydro- 
dynamic IF structure by Axford 1961). We shall see that, 
in this approximation, only fronts obeying the evolutionary 
conditions can have smooth structures. For the structures to 
be generic, the jumps across them must satisfy the strong 
evolutionary conditions, although singular classes of weakly 
evolutionary fronts with internal constraints on their flow 
structures are also possible. 

In this model the form of an IF is given by the variation 
of the roots of equation (jl^) with the temperature of the gas, 
specifled by a. The left hand side of this equation is a quartic 
independent of a, which is positive at 5 = 0, 00 and zero at 



a mcreasin 




0.5 



1.5 



Figure 2. Schematic plot of equation (^), where the solid curve 
corresponds to the l.h.s. while the dashed curves correspond to 
the r.h.s. for various values of a. The points where the solid curve 
passes through zero give the solutions of the isothermal shock 
jump conditions. 



5 = 1, and has two or four positive roots (since the condition 
that its value is zero is the normal MHD shock condition for 



isothermal gas, cf. Anderson 1963). The right hand side is 
a quadratic, which is zero at ni^S = 2ri and depends on 
Q through its (negative) scale. If the value of a increases 
steadily through the front, the manner in which the solutions 
vary can be followed, as illustrated by the schematic plot, 
Fig. I 

For the upstream conditions (a = 1), there will be either 
two or four solutions where the quartic curve crosses the 
axis. Of these, one is the trivial solution, 5 = 1, and at most 
one corresponds to an evolutionary shock. 

When internal heating occurs in an IF, the solutions 
at a given a correspond to where the solid curve in Fig. ^ 
crosses the corresponding dashed curve, —{a — l){nilS — 
2771). If gas enters the IF at the point marked o, in region 2 
between the Alfven speed and the fast mode speed, then it 
can move following the arrow as a increases. For sufficiently 
large a, the dashed curve becomes tangent to the solid curve: 
when this occurs, the gas is moving at the slow- or fast-mode 
speed, and the IF is called a critical front. In between the 
initial and final solutions, smooth, steady IF solutions must 
remain between the same characteristic speeds as they were 
when they started. The strong evolutionary conditions give 
just those cases in which a front structure calculated for a 
smoothly varying, monotonically increasing a (and non-zero 
perpendicular magnetic fleld) has a continuous solution from 
the upstream to the downstream case. 

When the upstream flow is at the fast or the slow critical 
speed, the l.h.s. of equation ( p^ has a second root 5 = 1. 
Thus for a > 1, this pair of roots disappears, and a forbidden 
region is generated. For a < 1 there is no forbidden region. 

For any a, the points where the solid curve crosses the 
dashed curve in Fig. ^ are related to each other by the 
isothermal MHD shock jump conditions, so a steady shock 
can form anywhere within the IF structure for identical up- 
stream and downstream conditions. However, since such a 
flow is over-specifled, the equality of the speed between the 
shock and IF is a coincidence which will be broken by any 
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perturbation of the flow (in which case the shock wiU escape 
from one or other side of the IF) . This situation is in direct 
analogy with strong R-type IFs in unmagnetized flows. 

Strong D-type fronts, for which the exhaust leaves the 
front rather above the critical speed, can occur where the 
heating is not monotonic. For these, the highest tempera- 
ture is attained when the flow is at the critical speed (i.e. 
the curves in Figure ^ become tangent exactly at the high- 
est value of a), and as it subsequently cools the solution 
can move back up the other branch. These fronts will form 
a more restrictive limiting envelope on the allowed weak 
solutions than that given by the critical solutions. The evo- 
lutionary conditions are necessary but not sufBcient, so this 
behaviour would be expected when more detailed physics 
was included. The actual envelope will correspond to the 
case for critical fronts with exhausts at the highest tem- 
perature attained within the front. Since unmagnetized IF 
models suggest that any overshoot in the temperature of the 
gas is likely to be small, the envelope will probably not differ 
greatly from that found for critical solutions, although the 
transonic nature of these fronts can be important in deter- 
mining the structure of global flows. 

Equation (|l^) suggests that no MHD flow can pass 
through the Alfven velocity (where 5ml = 2771) in a front 
unless it has zero perpendicular magnetic field. Heating the 
gas will generally move the fiow in regions 2 and 3 away from 
the Alfven speed in any case (see Figures ^ and ^ . However, 
where the perpendicular magnetic field is zero, a smooth, 
weakly evolutionary, front structure can be found (the in- 
ternal constraint being zero perpendicular magnetic field). 
The internal structure will be identical to an unmagnetized 
IF. Indeed, the strong-D hydrodynamic front can become, by 
analogy, an 'extra-strong' front which passes through both 
the Alfven and sound speeds (for zero perpendicular field, 
the slow and fast velocities are each equal to one of these). 

By analogy with equation (^o|), the y-components of 
velocity and magnetic field are zero everywhere if the MHD 
equations apply throughout the front (except if it were to 
pass through the Alfven velocity). Components in these di- 
rections can be generated if the velocity coupling between 
different components of the fiuid - electrons, ions , neu- 
trals or dust - is not perfect (as in shock structures, Pilipp 



rapidly away from the source, the speed of the front de- 
creases so that eventually the ionized gas exhaust is at the 
fast- mode speed (at the fast-R-critical velocity). When this 
occurs, two roots of equation ( pT| ) merge, and become com- 
plex for smaller mi. As a result, the front will then have to 
under go a transit ion of some sort. As in the unmagnetized 
case (Kahn 1954), if the size of the ionized region is large 



& Hart:iuist 1994; Wardle 1998). These components must, 
however, damp at large enough scales: far beyond the front 
the magnetic field must be in the same plane and of the same 
sign as the upstream field, from the evolutionary conditions. 
Exactly this behaviour has been found to occur in time- 
dependent multifiuid models of MHD shocks (Falle, private 
communication). A full treatment of ionization fronts in 
multicomponent material is, however, beyond the scope of 
the present paper. 



5 TIME DEVELOPMENT 

In this section, we discuss the development of the IF in a 
magnetized Hii region, by combining the well-understood 
development of IF in unmagnetized environments with the 
classes of physical roots of equation (|l^) found above. 

An IF driven into finite density gas from a source which 
turns on instantaneously will start at a velocity greater 
than the fast-R-critical velocity. Unless the density decreases 



compared to the lengthscales which characterise the inter- 
nal structures of shocks and IFs, the IF will evolve through 
emitting (one or more) shocks. The evolution of an initial IF 
discontinuity can be treated as a modified Riemann prob- 
lem because the speed of the IF is determined by the fiow 
properties on either side of it, together with the incident ion- 
izing fiux which we assume varies slowly. The development 
of this modified Riemann problem will be self-similar, just 
as for conventional Riemann problems. One complication is 
that the IF may be located within a rarefaction wave, but 
this does not occur for the circumstances we discuss in the 
present section. 

The simplest possibility for a slowing fast-R-critical IF 
is that it will become fast-D-type by emitting a single fast- 
mode shock. This has obvious limits to the cases where the 
magnetization is zero (where the shock is a normal hydro- 
dynamic shock), and where it is purely perpendicular. If the 
speed of the IF is specified by the mass flux through it, then 
the leading shock driven into the surrounding neutral gas 
must be a fast-mode shock, since the upstream neutral gas 
must still be advected into the combined structure at more 
than the Alfven speed after the transition, and so only a 
fast-mode shock can escape. 

While there may be no fast-D-type solutions at the 
value of ^1 which applied for the fast-R-type front, the fast- 
mode shock moving ahead will act to increase the value of 
^ upstream of the IF, since the (squared) increase in the 
perpendicular component of the magnetic fleld dominates 
over the increase in the gas density after the shock. At fast- 
R-criticality, there is in general a second solution with the 
same upstream and downstream states in which a fast-mode 
shock leads a fast-D-critical IF, because the l.h.s. of equa- 
tion ( p^ ) (for which a zero value implies an isothermal shock 
solution) is positive for large 5 and for post shock flow at 
the Alfven speed (i.e. mfd = 2r;i), unless the flow into the 
front is at the Alfven velocity, or rji or ^1 is zero. For ex- 
ample, for a downstream state 772 = = 5xl0~'^, there are 
two evolutionary solutions: a 1 — > 1 front with rji = 0.994, 
^1 = 0.248, mi = 20.0 and S = 0.503, and a 2 ^ 2 front 
with rji = 3.11xl0"^ ^i = 11.1, mi = 0.625 and 5 = 16.1. 
Figure |l| parts (c) and (b), respectively, contain points which 
correspond to these solutions. These two upstream states are 
linked by a fast-mode shock with the same velocity as the IF. 
Thus the emission of a single fast-mode shock is a valid so- 
lution of the modified Riemann problem which occurs as the 
fiow passes through criticality whatever the internal struc- 
ture of the shock and IF, so long as the evolutionary shocks 
and IFs exist. (Since the l.h.s. of equation ( p^ is greater 
than zero so long as r;i 7^ for 5 = 0, an equivalent argu- 
ment holds for slow-critical transitions.) 

It is possible that further waves may be emitted at the 
transition, for instance a slow-mode shock into the neutral 
gas together with a slow-mode rarefaction into the ionized 
gas. For the model resolved IF, this seems unlikely to occur 
unless the fiow velocity reaches the slow-mode speed some- 
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Figure 3. Plot of the dilution factor, S, within magnetized IFs 
for r]i = 3x10"^ and = 11 (cf. Figure ^(b) and the numerical 
example in the text), and decreasing towards 2r7i (with values 
for mi , plotted from top to bottom of the fan entering at <5 = 1, of 
mi = 1, 0.7, 0.5, 0.4, 0.35, 0.3, 0.27, 0.26, 0.252, 0.247, the critical 
solution and 0.22). Also shown dashed are the non-evolutionary 
Alfven solution for critical mi. The solutions terminate for a < 
100 for both high and low values of mi . Note also the development 
of the gradient discontinuity at a = ^i + l in the density structure. 



where within the fast-D-critical IF. These other solutions 
are not seen in our numerical simulations below. Additional 
solutions would also make the development of the IF non- 
unique, if the simpler possibility is allowed. 

The fast-critical transition may be followed using the 
jump conditions for a front with its exhaust at the fast-mode 
speed (i.e. a fast-critical front). In the smoothly- varying a 
model of Section ^, the upstream and downstream states 
can be joined by a front in which an isothermal MHD shock 
is at rest in the IF frame anywhere within the IF struc- 
ture, since the quantities conserved through the front are 
also conserved by the shock. The evolution of an IF through 
criticality will occur by an infinitesimally-weak fast-mode 
wave at the exhaust of the IF moving forward through its 
structure and strengthening until it eventually escapes into 
the neutral gas as an in dependent shock (as illu strated for 
recombination fronts by Williams & Dyson 1996). 

The escaping fast-mode shock leads to a near- 
perpendicular field configuration upstream of the IF. This 
boost in the perpendicular component is required if the tran- 
sition is toproceed though the fast-D type solutions, which, 
as Figure pj illustrates, are near-perpendicular (large ^i) ex- 
cept where is very large or very small. This will result in a 
rapid change in downstream parameters across a front where 
the upstream field is nearly parallel to the IF. As an exam- 
ple, the flow downstream of the IF will either converge onto 
or diverge from lines where the upstream magnetic field is 
parallel to the ionization front, and as a result may produce 
inhomogeneities in the structure of Hii regions on various 
scales (from bipolarity to clumps). 

As the velocity of the front decreases further, eventually 
it will approach the Alfven speed. We find that the evolution 
depends qualitatively on whether a — 1 is smaller than ^i. 
In Fig. ^, we plot 5 as a function of a for a range of values 
of mi (i.e. the internal structure of the fronts in the simple 
model of Section We take r]i = 3xl0~^ and ^i = 11, 
corresponding to our numerical example above, so the ratio 
of the upstream Alfven speed to the upstream isothermal 



sound speed is 0.245. For mi slightly larger than 0.245 (i.e. 
just super- Alfven), the curves remain flat until a is close to 
^1-1-1 and then turn upwards when they reach this value. If 
the value of a in the fully ionized gas were less than ^i -I- 1, 
the solutions would move through a case where the flow is of 
uniform density and moves at the Alfven speed throughout 
the front before becoming slow-R-critical for some ml < 2ri, 
as can be seen in the leftmost part of Fig. ^ for fronts in 
which the maximum a is smaller than 12. 

For larger values of a, however, the solutions develop 
a gradient discontinuity in their structure when the inflow 
is at the Alfven speed, mf — 2rj. At this discontinuity, the 
transverse component of the magnetic fleld becomes zero 
(i.e. switch-off occurs). Once the propagation speed of the 
IF drops below the Alfven speed, there is no form of steady 
evolutionary structure with a single wave in addition to the 
IF which is contiimous with that which applied before. Non- 
evolutionary 3^2 type solutions do exist for fronts just 
below this limit, but in numerical simulations (see Section^ 
these break up. A slow-mode switch-off shock moves into 
the neutral gas and a slow-mode switch-on rarefaction is 
advected away into the ionized gas. Between them, these 
waves remove the parallel component of magnetic fleld at the 
D-type IF between them. A precursor for the rarefaction is 
apparent in internal structure of the critical front, Figure ^ 
We find in numerical simulations that the IF which remains 
is trans- Alfvenic. Note that a steady resolved structure is 
possible for such (weakly evolutionary) fronts only because 
it has exactly zero parallel field throughout. 

Analogous processes must occur in an accelerating fast- 
mode IF with weak parallel magnetization when rji is very 
large: in Fig. ^ (d), it is the fast-D-critical line rather than 
the slow-R-critical line which meets the Alfven locus at finite 
6- 



6 NUMERICAL SOLUTIONS 

In this section, we present some numerical examples to illus- 
trate the processes discussed in the preceding sections. We 
have implemented linea r sch eme A as described by Falle, 
Komissarov & Joarder (1998) in one dimension, and added 
an extra conserved variable corresponding to the flow ion- 
ization. To study the local development of the IFs we have 
neglected recombination terms and just chosen to set the 
mass flux through the ionization flux as a function of time. 

Figure ^(a) shows the propagation of an IF into gas 
with density 1, = O.Sa/Itt, Bx = 2\/47r. We reset the flow 
temperature at the end of each time step so that p = (0.1-1- 
0.9a;)p (the temperature ratio between ionized and neutral 
gas is rather small so that any shells of shocked neutral gas 
are more easily resolved) . For these values, the characteristic 
speeds in the neutral gas are Ds = 0.046, Vs, ~ 0.3 and Vf — 2. 
In the first figure, the incident flux varies as 80/(20 -|-t), and 
the transitions through fast-R, fast-D and slow-D are clearly 
visible, while the slow-R stage (which is often narrow in the 
parameter space of Fig. ^ is less clear. 

In Figure ^b) the upstream conditions are the same, 
but the flux was set to a constant value of 0.26 in order to 
isolate the slow-R transition. This is between the upstream 
slow-mode and Alfven speeds, and equation (|l^) predicts 
that a slow-R transition exists which will increase the flow 
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Figure 4. Evolution of an oblique magnetized IF for (a) gradu- 
ally decreasing incident radiation intensity; (b) constant intensity 
chosen to show slow-R-type IF. In (a), the frames show log base 
10 of density (solid) and ionization fraction (dotted), offset by a 
constant between each plot, plotted at times 4, 8, 16, 1024. The 
topmost plot shows a fast-R IF, while third, fourth and fifth show 
a fast-D IF together with a fast shock moving to the right. In the 
next two, the IF changes to a slow-R type while in the final two 
a slow shock has moved off leaving a slow-D type behind. In (b), 
a more clearly resolved slow-R type front is shown. In this figure, 
the density is plotted linearly, scaled by a factor 0.1, and the plots 
are at times 200, 400, ...,2000. 



density from 1 to 1.18 (for comparison, the 3^4 jump rela- 
tions require the density to increase to 1.75). In the simula- 
tion, a weak fast-mode wave propagates off from the IF ini- 
tially, but does not greatly change the upstream conditions. 
The slow-R IF which follows it increases the density, and is 
backed by a rarefaction because of the reflective boundary 
condition applied at the left of the grid. The small over- 
shoot within the front is presumably due to numerical vis- 
cosity, and can be remov ed by broadeni ng the IF (e.g. by 
the method described in Williams 1999). The plateau be- 



tween the rarefaction and the IF has density 1.22, which is 
in adequate agreement with the jump conditions (particu- 
larly when account is taken of the slight perturbation of the 
conditions upstream of the front). 

These numerical solutions illustrate the orderly progres- 
sion of a magnetized IF through the various transitions de- 




Figure 5. Evolution of an IF from an initially non-evolutionary 
solution. The curves show the density (0.1 times the log base 10, 
solid), ionization fraction (multiplied by 0.1, dotted) and 0.1 X 
Bx/V^ (dashed), plotted at times 400, 800, 1200, 4000 with 
an offset between each plot. 



scribed in the previous section. In a realistic model, however, 
the density perturbations generated by the transitions will 
have important effects on the evolution, as the consequent 
changes in recombination rates alter the flux incident on the 
IF. These processes should ideally be studied in the context 
of a two- or three-dimensional global model for the evolution 
of magnetized Hii regions, which is beyond the scope of the 
present paper. 

In Figure ^, we illustrate the development of the flow 
from initial conditions in which a non-evolutionary IF is 
stationary in the grid. The upstream (neutral) gas has den- 
sity 1, Bz = 0.3-\/47r, and Bx = \/47r and moves into the 
front at Vz = —0.253 (with no transverse velocity), while the 
downstream (fully ionized) gas has p = 0.423, Vz = —0.598, 
vx = 1.688 Bz = 0.3^4^, and B^ = -0.424 so the IF is 
of 3 — > 2 type. We set the pressure as above, and the value 
of the incident flux as 0.253 so the initial IF would remain 
steady in the grid. The IF breaks up immediately, driving 
a slow- mode shock away to the right, into the neutral gas, 
while a slow-mode rarefaction moves away to the left, into 
the ionized gas. The D-type IF which remains is marginally 
trans- Alfvenic, but has zero transverse magnetic fleld. 

To study this further, in Fig. |^, we illustrate a simu- 
lation of an IF close to the (hydrodynamic) D-critical con- 
dition, with a parallel magnetic field which makes it trans- 
Alfvenic. When perturbed with small but significant per- 
pendicular field components, this IF again switches off these 
fields by emitting slow-mode waves. It remains stable, sat- 
isfying the hydrodynamic jump conditions as a weakly evo- 
lutionary solution of the MHD jump conditions. 



7 CONCLUSIONS 

We have presented the jump conditions for obliquely- 
magnetized ionization fronts. We have determined the re- 
gions of parameter space in which physical IF solutions oc- 
cur, and have discussed the nature of the interconversions 
between the types of front. Fast-D and slow-R solutions with 
high transverse fields are found in the region of front veloci- 
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Figure 6. Evolution of an IF obeying tlie liydrodynamic D- 
critical conditions at rest, perturbed by an initial magnetic field. 
The upstream density is 1 and velocity is —0.0513, while the 
downstream density is 0.0513 and velocity is —1. A parallel field 
of Bz = 0.1\/47r means this front is trans-Alfvenic. The IF 
is perturbed with a perpendicular magnetic field of magnitude 
0.03\/47r, in different planes up- and down-stream. The solid line 
shows the density, the dotted line is the ionization fraction and 
the dashed line is one component of the perpendicular magnetic 
field (which has been scaled to a maximum amplitude of 0.5). 
Also shown are the positions of the characteristic waves from the 
initial discontinuity (F is fast, A is Alfven, S is slow, M is mass), 
based on the initial states. The waves which change the magnetic 
field are, from left to right, a slow-mode shock plus an Alfven 
wave, a slow-mode switch-on rarefaction (which switches off the 
field as it is being advected away from the IF) and a slow-mode 
switch-off shock propagating into the upstream gas. The front 
which remains is trans-Alfvenic, but has effectively zero perpen- 
dicular field. 



ties forbidden by the hydrodynamic jump conditions: in the 
evolution of an H ii region, the fast-mode shock sent into the 
neutral gas by the fast-critical transition will act to gener- 
ate these high transverse fields. In the obliquely-magnetized 
case, the fronts are significantly perturbed as long as the 
Alfven speed in the neutral gas is greater than cf /2c2 • How- 
ever, the stability of parallel-magnetized weakly evolution- 
ary IF means that the flow may still leave at the isothermal 
sound speed in the ionized gas over much of the surface of 
magnetized globules exposed to ionizing radiation fields. 

Large (<; 100 /iG) , highly ordered magnetic fields have 
been obser ved in the molecula r gas surrounding som e H ii re- 



gions (e^ 



Roberts et al. 1993 



Roberts et al. 1995). Roberts 



et al. (1995) suggested that the highest observed magneti- 



zations in S106 are associated with unshocked rather than 
shocked gas, as a result of the relatively low density and ve- 
locity shift observed for the strongly magnetized gas. They 
suggested that the magnetic field becomes tangled close to 
the IF, leading to the decrease in detectable magnetization. 

It is interesting to compare these results with the ex- 
ample fast-critical front we discuss in Section ^. Scaling the 
parameters of this front to an exhaust hydrogen density of 
10*cm~'^, typical of an ultracompact Hii region, a mean 
mass per hydrogen nucleus of 10"^* g, and a sound speed in 
the ionized gas of lOkms"^, the limiting 1^1 front takes 
the flow from an upstream state with nu = 5.03x10^ cm~^, 
V = (-0.0497, 0, 20.0) kms'^ and B = (17.7, 0, 35.4) ^iG to 



an exhaust with nu = 10* cm ^, v = (0,0, 10.1) kms"^ and 
B = (35.4, 0, 35.4) /iG (while the fast-R-critical speed is al- 
most exactly twice the exit speed of the front, this ratio 
decreases in more strongly magnetized fronts). The limit- 
ing 2 — » 2 front takes the flow from nn = 1.61x10^ cm~'^, 
V = (1.78, 0, 0.625) kms"^ and B = (671, 0, 35.4) /xG to the 
same flnal state. A fast-mode shock ahead of the front and 
at rest with respect to it would change the upstream state 
from that of the limiting 1^1 front to that of the limiting 
2~*2 front. 

This fast-mode shock, which precedes a limiting fast- 
weak-D type front, boosts the x-component of the magnetic 
fleld from 18 fiG to 670 fiG, with a 20kms~^ change in the 
2- velocity and only a 2 km s~^ change to the a;- velocity com- 
ponent. A succeeding slow-mode shock would further in- 
crease the z- velocity component and gas density while weak- 
ening the a;-component of magnetic fleld, without recourse 
to fleld-tangling. If we tentatively identify OH component 
B of Roberts et al. as fast-shocked material and component 
A as doubly-shocked material, component B is more edge- 
brightened and has a smaller blueshift than component A as 
would be expected. Component A is kinematically warmer 
and most blue shifted towards the centre of the region. The 
strong line-of-sight magnetic fields in S106 are seen at the 
edges of the region, in a 'toroidal' distribution. The (poorly 
resolved) line-of-sight velocity of component B has little gra- 
dient in the equatorial plane of the region, but this might 
result in part from a combination of fiow divergence and the 
value of _Bios (which is measured close to the centre of the 
region) being rather larger than that of in our example. 
While the qualitative properties of this assignment are at- 
tractive, it remains to calculate a proper model tuned to the 
properties of the region, in particular its geometry. Never- 
theless, the present discussion at least illustrates how both 
fast and slow shocks should be considered in the analysis of 
regions with well-ordered magnetization. 

In future work, we will model in detail the global struc- 
ture of magnetized Hii regions and the local structure of 
photoevaporated magnetized clumps. 
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